Small Black Holes on Branes: Is the horizon regular or singular ? 
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We investigate the following question: Consider a small mass, with e (the ratio of the Schwarzschild 
radius and the bulk curvature length) much smaller than 1, that is confined to the TeV brane in 
the Randall-Sundrum I scenario. Does it form a black hole with a regular horizon, or a naked 
singularity? The metric is expanded in e and the asymptotic form of the metric is given by the weak 
field approximation (linear in the mass). In first order of e we show that the iteration of the weak 
field solution, which includes only integer powers of the mass, leads to a solution that has a singular 
horizon. We find a solution with a regular horizon but its asymptotic expansion in the mass also 
contains half integer powers. 
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O ' I. INTRODUCTION 

o : 

<N . ..... n 

Black holes in theories with extra dimensions have been studied widely. Myers and Perry [1| found Schwarzschild 
type solutions (MPS) in D-dimensional asymptotically flat space. Black hole solutions were also found in asymptoti- 
cally AdS space 0, Q No non-trivial closed form black hole solutions, other than the black string solution jj|, which 
extends in a uniform manner from the brane into the extra dimension, have been found in three-brane theories of the 
Randall-Sundrum type Given that there is considerable interest surrounding the production of black holes at 

accelerators |7j , and in collisions of cosmic rays [8fl , it is important to develop approximate methods to find black hole 
solutions in Randall-Sundrum brane world theories. 

1^ , Some initial attempts at finding black hole solutions centered on deriving the induced metric on the brane by 
solving the Hamiltonian constraint conditions |9( . Some of the induced solutions do not arise from matter distributions 

Q ! confined to the brane 0] . Linearized solutions about RS backgrounds 0] as well as numerical solutions [HE! have 

Tj" ■ also been derived. 

In a recent paper [3] we have studied the metric around a small mass that is confined to the TeV brane in Randall- 
Sundrum type I scenario (RSI)[j|. The configuration is characterized by two length scales; the Schwarzschild radius 
fi, which is related to the mass /i = ^/8GsM/ (Sir), and the curvature length of the bulk £, which is related to the 
cosmological constant and the brane tension. We study the case \x <C £ where the metric can be expanded in the 
dimensionless parameter, e = fi/£ <C 1. There are two ways to make e < 1: (1) /i — > while £ is finite. This limit 
$H . is just the original Randall-Sundrum scenario in the absence of matter. Keeping £ finite and expanding in /i (or 
actually in M) is called the weak field approximation. First order in M is linearized gravity. Linearized gravity is 
, valid at distances much larger than the Schwarzschild radius, r 3> /i. (2) £ — > oo while fi is finite. This limit is a 
five dimensional black hole in an asymptotically flat background. The metric for this configuration is the Myers and 
Perry solution (MPS). We studied the expansion in e with /i kept finite. We call the first order expansion in e - the e 
solution. The e solution is valid at distances much smaller than the bulk curvature length r ^ £. 

On one hand, the e solution is needed to study the horizon and the thermodynamics of the black hole, for which the 
linearized solution is unsuitable. On the other hand, the e solution cannot be fixed uniquely without satisfying certain 
boundary conditions. We need the linearized solution to identify the mass and to satisfy the junction conditions on 
the Planck brane because the distance between the branes in RSI is of the same order as the bulk curvature length 

Since we need both the e solution and the linearized solution we look at the region where both solutions are valid, 
^ < r < 1 In this region we expand the linearized solution to first order in £~ 1 and the e solution to first order in 
M. We require that these two expansions will coincide. Using this method, we are able to incorporate the asymptotic 
characteristics into the short-ranged e solution. In other words, the linearized solution is the asymptotic boundary 
condition for the e solution. 

In this paper we study the regularity of the horizon. In section [D] we summarize the results from In section 
IIIII we analyze the horizon and formulate the conditions such that the e solution will be regular on the horizon. We 
also calculate the thermodynamics parameters of the black hole assuming it is regular. In section Hvl we represent the 
e solution with Legendre functions. We assume that in the asymptotic region the metric is given by the linearized 
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solution (linear in the mass). Another reasonable assumption is that when iterating the linearized solution, the post 
linearized solution includes only integer powers of the mass. This assumption is motivated by the post Newtonian 
behavior of 4 dimensional black holes, where the asymptotic solution is an expansion in integer powers of G^M/r. 
We show that under this assumption the horizon is singular. In section we show that one can construct a regular 
solution provided that the post linearized solution includes half integer powers of the mass. We find it interesting 
that the first law of black hole thermodynamics, dM = T dS, eliminates the terms of order M 3 / 2 from the asymptotic 
expansion. 



II. THE e SOLUTION 



We use the following ansatz for the metric 

A(p,ip) , , 



ds l 



i(p, ip) 



-B(p,ip)dt 2 



dp 2 + 2V(p, ip)dp dip + p 2 U(p, ip)dip 2 + p 2 sin 2 iP dfl 



(1) 



where dtt?, is the metric on the unit 2-sphere. The conformal factor is taken from RSI as a(p,ip) — (1 — ep\ cos-0|) -2 . 
The TeV brane is located at ip = n/2 as can be read off of the conformal factor. We work only in the interval 
< ip < 7r/2, and assume Z 2 symmetry about the brane. The point mass is located at p = 0. The e solution is just 
MPS solution with e corrections 



1 



1 



A(p,i>) = l + eAx(p,iP) , 
U(p,ip) = 1 + . 
V(p,il>) = eVrip,^) . 



(2a) 

(2b) 
(2c) 
(2d) 



We have rescaled the coordinates with p such that p — 1 corresponds to the MPS horizon at the Schwarzschild radius. 
The gauge is partially fixed by the coefficient of dfl 2 ,. There is still a gauge freedom associated with the coordinate 
transformation p — ► p(l — eF), ip — ► ip + eF tan ip. 

The e solution can be written in terms of the gauge function F(p, ip) and the wave function H(p, ip). The solution 
is E3 



Bi 
At 



P 

8p(p 2 - 1) cos i\> - 12H(p, ip) - 4 tanipH^p, ip) + 2F(p, ip) + 2pF p {p, ip) 
2 tan 2 ip 



P 



[6pH(p, ip) - (2p 2 - l)H p (p, ip) - pF(p, j>) - pcot ipF^p, ■ 



-2p 2 (2p 2 - l)sinV> + 2p 2 tanipH tP (p : ip) 



+ 



p(2p 2 - l)tan 2 ipH^jp, ip) 
p 2 - 1 



p 2 taxitpFJp, ip) 



p 3 F^(p, ip) 



The wave function H(p, ip) satisfies the differential equation 

1 



(p 2 - 1) ( H pp - -H p 



H 



ijjljj 



cos 2 ip + 1 

2 ~ — ; r-". 1 

sm ip cos ip 



(3a) 
(3b) 
(3c) 

(3d) 
(4) 



Junction conditions on the brane. Israel junction conditions 15] on the brane are simply Bi^(p,ir/2) = 0, 
Ai t ip(p,ir/2) = 0, Vi(p,ir/2) = 0, and Ui(p,ir/2) < oo. These conditions imply the following 



i^(p,7r/2) = , 

(W + tampH^), =-2p(p 2 -l), 

,v ip—n/2 

H^(p,n/2) = , 

F(p,ir/2) = 6H(p^/2) - ^—l H . p (p^/2) 

P 



(5a) 
(5b) 
(5c) 
(5d) 
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Asymptotic boundary conditions. As mentioned earlier, we take the linearized solution as an asymptotic boundary 
condition. As a result, the asymptotic form of the wave function is 

3 

H(p,ip) = (5 H (40- sin 40) + ^—(4 sin 20 - 80 cos 20 - 120 + 3 sin 4*0) 

16 sin 64 sin 

+ — (1024(5 + 80 + 32-0 cos 20 + 120cos40- 16sin20- 5sin40) + 0(p- (1+2a) ) , (6) 

1024p sin 

where <5 and So are arbitrary constants. The asymptotic form of the gauge function is 

F(p, 0) = t ( cos + J-) + 6S + + (p-^) , (7) 

2 \ sm0/ p 

where F 2 (tt/2) =85- 5tt/64 and Fg(n/2) = 47tt/64 - 245. 

In the functions 10 and J7J we have neglected corrections of order p~( 1+2a ) where a > 0. The reason is that such 
terms will contribute to order M 1+a , which is beyond the linearized solution. 

A reasonable assumption is that if one will go beyond the linearized solution, only integer powers of M will appear 
in the metric. This assumption will force the functions H(p,ij}) and F(p,tp) to be odd functions of p, therefore one 
must set 5q = and the neglected terms are of order p~ 3 . 



III. ANALYZING THE HORIZON 



The configuration is called a black hole if there exists a Killing horizon, i.e. a null surface with a Killing field normal 
to the surface. The ansatz Q is static, and the static Killing vector becomes null on the surface B = 0. The static 
Killing vector is £ M = 5$ , it is normalized such that asymptotically £ 2 = — 1. 

A vector normal to the surface B = Q must be proportional to = B >fi . If the surface B = is null, then the 
normal vector must be null on B = 0. 



Therefore, the derivative B j, must vanish on the surface B = 0. This implies that the surface B(p, 0) = is actually 
defined as p = pn = constant. The ansatz for B{p, 0) is 



BM)=(l-^jbM), (9) 

where b(p,ip) does not vanish at p = pn, but it is regular such that lhn p —> PH (fi — pH)b(p,ip) = 0. 
The surface gravity, k, for the ansatz is defined as 



= - ^V'^aW 



AA 



(10) 

B=0 



According to a theorem by Racz and Wald fl6| , if the surface gravity is not constant on the horizon, the horizon is 
singular. This can be verified in the ansatz Q by calculating the Kretchmann scalar around the horizon. Provided 
that the function B(p, ip) is given by Eq.@, the Kretchmann scalar can be expanded in powers of B 



B(p,ip) p 2 a 2 U 



R" vXa R^ = sr. „/.\ l:tr + re s ular terms ■ ( n ) 



Clearly, at any point where ^ 0, the surface B — is singular. If the surface gravity is not constant but vanishing 
at some points on the horizon, these points will be regular. An ansatz for the function A(p,ip), which is compatible 
with Eas.pHO)!. is 

A M = b ^l + ( 1 -^) aM) , (12) 

where lim p ^ Pff (p - p H )a(p, 0) = 0. 
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e expansion on the horizon. Take the ansatz Q together with Eas.(|9l H"2*) l and expand in e. Zero order in e should 
be MPS solution, therefore, the following e expansions should be used 

b(p,il>) = l + eb 1 M) , (13a) 

a(p,ip) = eai(p,tp) , (13b) 

p H = l + <, (13c) 

K=l + e X . (13d) 



Comparing Eas. (|13|l with Eqs.J2J), one can solve for b\ and a 



2(~p 2 B 1 (p,j;) 

h{p,V) = 5 — i ' ( 14 ) 

p z — 1 



P A — 1 



The expansion in e on the horizon is possible only if the functions b\ and a\ are finite at p = 1. Therefore, the 
numerators in Eas. i|14l I15fl must vanish at p = 1, and one can calculate Q, 6i(l, V0j an d X 

C = limiBi(p,VO , (16) 

hM) = -lin f I "< P , (17) 

X = limLli(p,VO-&i(l.VO-C. (18) 

Both £ and x are constants, therefore Eas. 1)161 118[) put restrictions on the solution ©. Equation (|16fl implies that 

C = -F(l, %l>) = constant F^(l, tp) = . (19) 

Substitute the solution in Eq. l|18|l to calculate the correction to the surface gravity 

X = 2(3H(l,ip)+taD.iiHrf(l,ip)) = constant . (20) 

So, the combination 3H(1,iJj) + tanip H ^ (1, ip) must be constant. 

Another constraint comes from Ea. (|3d|l . The function V\ must be finite at p — 1 (just like b\ and a±). This implies 
that 

lim(p 2 -l)Fi(p,V) =tan 2 V^(l,V')+^(l,^) = . (21) 
p->i 

Since F^(l,-0) = we deduce that H^(l,ip) — 0. The functions and -F are constants on the horizon and using 
Ea. (|5d(l can be evaluated as 

H(l, if,) = H(l, tt/2) , F(l, VO = F(l, tt/2) = 6ff (1, tt/2) - p (l, tt/2) . (22) 

To summarize, if the e expansion is valid on the horizon and the horizon is regular then 

i^(l,V) = 0, (23a) 

^(1,V) = 0, (23b) 

= 1 + e( = 1 + e[-6H(l, tt/2) + ff p (l, tt/2)] , (23c) 

«= l + ex= l + 6eff(l,7r/2) . (23d) 

Thermodynamics. If the horizon is regular then one can talk about thermodynamics of the black hole. The zeroth 

law of black hole thermodynamics states that the surface gravity is constant on the horizon The temperature, 
T, associated with the black hole is 

1 = ^ = 27^(1-6^(1, tt/2)) . (24) 

1 K 
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One should remember that the physical dimension of the surface gravity is length 1 therefore it is rescaled with p, 1 . 
The entropy of the black hole is related to the area of the horizon 



S = t^- = 7^- / dip [ dff [ d<j> y/g^geeg^ 
4(_z 5 4(_r 5 j j a j 



27TP fPH 



3„3 /•tt/2 



2 1 + 6/2^(1,^) 

a-0 sm ip - 



p=1+eC G 5 Jo "" r r (1-ecos^) 3 



(25) 



Using the solution © and the constraints l|22l I23[l . one can verify that on the horizon Ui(l,ip) 
2tan 2 ?/> [£T p (l,7r/2) - H tP (l,ip)]. The entropy is 



S = 



2ttV 









hi 













tt/2 



(26) 



The integral in Eq. (|2(j[l can be simplified using integration by parts as 



tt/2 



dip sin 2 ip tan 2 ip [H iP (l,ir/2) - H tP {l, ip)] 



71/2 sin?/' . 3 



dip 



cos 2 ip 



sm 6 iP[H p (l,n/2)-H >p (l,iP)] 



sin ip 
cos ip 



{#,„(!, tt/2) - 



3tt 



t/2 r7r / 2 



tt/2) -— # p (1, tt/2) + 



dip 

o jo cos^ 
/2 # 



{sin 3 i\> [H p (l, tt/2) -ff p (l,^)]} 







cos 



[sin 3 ^,(1,0)] „ , 



(27) 



where in the last step we have used L'Hopital's rule to evaluate the boundary term at ip = tt/2. Using Equations 127|) 
and l|23c|) . the entropy (|26[1 can be written as 



2ir 2 p 3 
4G 5 



S = 



1 + ei -18i?(l,7r/2) + - 

7T 



l + #^(l,7r/2) + 



tt/2 



COS ?/> 



sm 



i>H p (i,ip)] 



(28) 



The first law of black hole thermodynamics states that T _1 = dS/dM. In all known black hole solutions the first 
law is satisfied as a result of Einstein equations. In the e solution it is not trivi ally satisfied , but should be imposed 
as a boundary condition. The mass appears in the entropy only through p, = y8GsM/ (37?) and e = p/£, so 



1 dS „ 



l + eJ-24ff(l,7r/2) + ^ 



,7r/2 dii 

l + F p ^(l,7r/2)+ / -5- [sin 3 Vff,p(l,'. 



/q COS'0 

Equations l|29|l and (|24|l must be consistent, therefore, there is another constraint on the function H 



(29) 



H(l,n/2) = 



27-k 



/ 2 dih 

1 + ^(1,71/2)+/ — ^ [sin 3 ^ff,„(l,', 



(30) 



IV. THE SINGULAR HORIZON 



A detailed study of the horizon in the e solution requires a specific representation of the function H(p,ip), which 
solves the differential equation 10} and the boundary conditions 10 H3 123|) . 

A possible representation for the solution is a combination of Associated Legendre functions in p and Hypergeometric 
functions in ip 



(31a) 
(31b) 
(31c) 



H(p,iP) = j dAi?(p;A)*(0;A) 

R(p;X) = p^p 2 - 1 \a{X)Q\^ 1 (2p 2 - 1) + b(\)PL± (2p 2 - 1) 
*(V>;A) = c(A) 2 F!( 



1-A 1 + A 5 . 2 



d(X) 



2 '2' 



sin -0 



-2-A -2 + A 1. 2 ,, 
'^^'~2' Sm ^ 
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A lengthy discussion about the expansion 1|31[) appears in |T3 | . There we assumed that the asymptotic solution beyond 
the linearized solution includes only integer powers of M. Therefore, the function H is odd in p and is given by 



H(j>,i/>) = py/p 2 - 1 



2 ! 2 1 , 3[sin(2^)-2^cos(2^)] 

— Q_ 1/2 (2p -1) 



+Q} /2 (2p 2 -l) 



4 W-l) 



1 3[Aip - sin(4t/j)] 
3tt 32 sin 3 V 

2 9Qi_ 1/2 (2p 2 - 1) 



sin 

d 2 7T 3 [8^ + 40 cos(40) - 3 sin(40)] 
sin 3 ^ ~ 3 



oo 

+ ^a(A = 2n)Q„_ 1/ 

n=2 

The set a(A = 2n) n > 2 is fixed by the junction condition i|5b 

a(A = 2n)„> 2 = 



32tt 2 sin 3 ip 

3[40> - sin(4V0] 
32 sin 3 ip 



l fr, 2 3 [71 cos(2n.0) sin(20) — cos(2ip) sin(2n^)] 

n-l/2(2/' — 1) — — ; 3 

8n(l — n^j sin ip 



16 



3tt(A 2 - 4) 



The parameter d 2 is undetermined. 

The problem appears when one tries to apply the conditions (|2UI I23fl to the solution (|32|l . 

• The associated Legendre functions of the second kind are not analytic at p = 1 

1 



(32) 



(33) 



pVp^IQI^P 2 - 1) ~ - 5 - + !)(P 2 - 1) In IP 2 - !| + <V - 1) • 
As a result, one can evaluate the function H(l, ip) but not the derivative H „(1, ^). 

The surface gravity, (|2U[1 . is not constant. The combination 3_ff + tan ipH j, is not constant on the horizon 



(34) 



K =l + 2e [3H(l,ip) + tan </-(!»] = 1 + 2e 



20 sin 3cos"0 3sin(2n0) 
1 } a(\ = 2n)— — : — — 



2n 



n=2 



= 1 + 2e smip . 

(35) 

As was mentioned earlier, according to the theorem by Racz and Wald [l6f . if the surface gravity is not constant on 
the horizon, the horizon is singular. One may notice that on the brane K,*/>l^ =ir /2 = 0. So, according to Ea. (|ll|l . 
although the horizon is singular, the singularity is naked from the bulk but it is covered on the brane. 



V. GIVE UP THE POST LINEARIZED ASSUMPTION 



We give up the assumption that the asymptotic solution can be expanded in integer powers of M, but still keep 
the linearized solution as the asymptotic boundary condition. According to Eq.|@ we can include terms of order 
It is convenient to work with the Legendre expansion (|31() . The asymptotic behavior of the associated 



P 



-(l+2a) 



Legendre functions of the first (second) kind is p\/ p 2 — lP(Q)^ x _ 1 y 2 (2p 2 — 1) ~ p 1+ ( ) A . Therefore, we can include 
associated Legendre functions of the second kind with A > 2 in the expansion. Actually, we have already included 
even A terms in the expansion (|32|l . We first try to include also odd A terms. 



H(p,i/>) = S + p^Jp 2 - 1 



— O 1 f2j 2 , 3^(20)- 20 cos(20)] 
3tt 1/21 P ' 8sin 3 



-1^(2^-1) 



5>(A)Q[ A _ 1)/2 (2p 2 -l 

A=3 



3tt 32 sin 3 ip 

2 dQ\_ 1/2 {2p 2 -l) 



d 2 tt 3[80 + 4-0 cos(40) - 3 sin(40)] 
sin 3 ip 3 



On 



32tt 2 sin 3 
3[40 - sin(40)] 
32 sin 3 ip 



3[A cos(A0) sin(20) - 2 cos(20) sin(A0)] 
2 A (4 - A 2 ) sin 3 



(36) 
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The free parameters are d 2 , So, and a(A = 2k + l)fc>i- Our goal is to make the surface gravity constant at p = 1. 
Adding the odd A terms does not change the junction condition (|5b|) and therefore the coefficients a(A — 2n) are the 
same as in Ea . (133(1 . We will try to fix the coefficients <5o, a (A = 2k + 1) such that the combination 3H + t&nipH^, 
which appears in the surface gravity (|2U[) . is constant. From Ea. l|36l) we find that 



3H(1, 0) + tan0# v ,(l,0>) = 3<5 + sin0 - ^ a(A = 2k + = ;-;//, L , T /-_>) 



fc=i 



2(2fc + l)sin0 



(37) 



To solve for a(2k + 1), one can multiply l|37[l by sin0 and expand in the set {s'm(2k + l)0} fc=o 1 , which is complete 
and orthogonal over the interval [0,7r/2] provided that the boundary conditions are /(0) = and /'(7r/2) = 0. The 
coefficients are 



6 a = H(l,ir/2) 
a(A = 2k + l) fc ; 



9tt 



-16 



(38) 
(39) 



3tt(A 2 - 4) ' 

Since we have multiplied l|37|l with sin0 one should verify that (|37|l holds for "0 = 0, and indeed it is. 

At this stage we have to verify the consistency of; (i) the representation of the function H , Eq. (|36(l , (ii) the zeroth 
law of thermodynamics, Eas. ((38l I39[l . and (iii) the first law of thermodynamics, Eq. (J3DI ■ 

We use the following expansions for the Legendre functions around p = 1 



^yV^LP.W - 1) = v(y + l)(p 2 1) 



(p 2 -i) 2 + eV-rr 



2 2 
^> + l) 2 



(p 2 - 1) [1 - ln(p 2 - 1) - 2 7E - 20(^ + f ) 
hv 2 + 5v + 2 



2v(v + 1) 



- \n(p z - 1) - 2 lE - 20(^ + 1) 



(40) 



+ 0(p 2 -lf ,(41) 



where is Euler's constant and ip(z) is the logarithmic derivative of the Gamma function, ip(z) = d/ dz\nT{z). 
First, we evaluate Ea. (|36jl at p = 1, tp = tt/2 



H(l,n/2) 



. 1 1 d 2 3 1^2, 



= ^0 + 



32 

rf 2 



2 32 

E 



(-1) A 16 3sin(A7r/2) 
2 ^ 3^(A 2 -4) A(4-A 2 ) 



(-l) fe 



tt^ (2fc + l)((2fc+l) 2 -4) 2 9tt 2 



(42) 



Comparing Ea.l)42[l with the result of the first law, Ea. l|38(l . we conclude that d 2 = 0. 

Second, we would like to evaluate the first law, Ea. (|30|l . For this we need to evaluate £f p (l,0). We use the 
expansions (gUl ED in Eq- GSJ and find that 



#p(1,0) 



sin 3 



[\n{p 2 - 1) + 2 1E + 20(1/2)] + hxty) . 



-70 + 40 cos 20 + 30 cos 4-0 — 2 sin 20 + sin 40 



32tt 



2 ^ (-1) A (A 2 - 1) 



E 



[A cos A0 sin 20 — 2 cos 20 sin A0] . 



50 — 2 sin 40 + 30 cos 40 



4z 



8ir sin 
(-1) A (A 2 - 1) 



*ZZ A ( A2 - 4 ) 2 







A + l 



A cos A0 sin 20 — 2 cos 20 sin A0 
sin 



(43) 



(44) 



(45) 



Let us deal with the function ho(ip), which multiplies the diverging term ln(/? 2 — 1). One can verify that ho(tp — 
7r/2) = and ho(ip = 0) = 0. Therefore, one can expand the function ho(tp) in the set {sin2n0} n=1 2 , which is 
complete and orthogonal over the interval [0, 7r/2]. It is easy to verify that for n > 1 



r/2 



ho (0) sin 2n0 dip = 



(46) 
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So, the function ho(ip) = and the first derivative of the function H is finite on the horizon. The first law, Ea. (|30|l . 
together with Ea. l|38|) can be used to find the value of <5 



60 = 



27-k 



h'^/2) 



zr/2 



dip 
cos 4> 



[/ii^sin 3 if\ 



From Ea. (|45|l one can deduce the following 



[/irWsin 3 ^] = 



5-5cos#- 12V>sin4V> 4 ^ (-1) A (A 2 - 1) 



-I 2 di; 



A(A 2 -4) 



A+l 



j/)=7r/2 

5 2 



. 



^ [/nW sin 3 ^]. = - + -E 7^77 $ ( n + ^ - ^ (V2)] = 2 • 
cost/; l j 'V 7r 7r z ^ (n^ — 1) 



(47) 

sin 2tp sin A^ . (48) 
(49) 
(50) 



Now, substitute in Eq.JJJJ to find Sq = 0. 

To summarize, if we give up the post linearized assumption and allow for asymptotic behavior of M n l 2 then we can 
satisfy the requirement that the surface gravity is constant and the solution appears to be regular on the horizon. The 



function H and the first derivative H p are finite on the horizon. 



The derivatives and H ^ are vanishing there. 



Therefore, based on Eq.|0J we can conclude that H pp is finite on the horizon as well. This is a sufficient condition for 
the metric to be G 1 . In general this is not sufficient to ensure that the the metric is not singular. To double check, 
we look at the Kretchmann scalar up to first order in e 



K = R abcd R a , 



bed 



72 

78 



576 

„8 



( F(p, ip) + p{p 2 - 1) cos + 



96 f 3H(p,iP) 
P \ 



taatpH^(p, ip) 



(51) 



As one can see, the Kretchmann scalar depends only on the first derivative H iP (p,ip), and therefore it is finite. An 
ambiguity is left for the second order in e. The Kretchmann scalar at the second order in e depends on higher 
derivatives of H, which might be divergent. However, second order in e in the Kretchmann scalar depends on the first 
and second orders in the metric. So nothing can be done at this stage before finding the second order corrections to 
the metric. These should cancel the divergence part coming from the first order in the metric. We should emphasize 
that the situation with the surface gravity is different. It is true that if the surface gravity is not constant in first order 
in e, it will affect the Kretchmann scalar only at second order in e. However, this contribution cannot be compensated 
by second order corrections to the metric. Therefore, the surface gravity must be constant at all orders. 

What about the asymptotic behavior of the metric? The inclusion of half integer powers of the mass changes the 
post Newtonian potential. For large p, the terms with a(A = 2k+ l)fe>i contribute to the function H (and throughout 
Ea. (|5dJ) to the function F) terms of order p~ 2k . As a result, the gravitational potential, g tt , acquires terms of the 
form 



a(A = 2fc + l) a(A = 2k + 1)(G 5 M) 



fe+3/2 



„2fe+2 



p2k+2 £ 



(52) 



where p — pp is a dimensionful coordinate. The term with $o would have contributed a term of the form 
S (G 5 M) 3 / 2 p~ 2 £^ 1 . But, since S Q = the half integer terms start only from M 5 / 2 . 

We learn that, in general, the asymptotic solution should be expanded in powers of (G5M) 1 / 2 . The lowest order 
must be (G5M) 1 and not (G5M) 1 / 2 otherwise there will not be a well defined conserved mass. One might expect that 
by iterating the lowest order there will be only integer powers of M. However, the boundary conditions on the horizon 
require also half integer powers of the mass. In particular it is the zeroth law of black hole thermodynamics (constant 
surface gravity) that forced us to include the half integer terms. The coefficients of the terms (G^M) k / 2 , k — 5, 7 . . ., 
which are a(A = 2k + l)fc>i, are completely determined by the zeroth law (|39|l . The coefficient of (G5M) 3 / 2 , which 
is 5q, is not fixed by the requirement that the surface gravity is constant. The actual value of the surface gravity 
depends on Sq. The first law of black hole thermodynamics, dM — T dS, forces So to vanish, and therefore eliminates 
the term of order (G5M) 3 / 2 from the asymptotic expansion. 

The dominant part of the post Newtonian potential remains (G5M) 2 /(r 2 £ 2 ), which is responsible for the precession 
of perihelion calculations. Since we don't know of any experimental evidence for higher order terms in the potential, 
it looks like the half integer terms cannot be detected experimentally. 
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VI. SUMMARY AND DISCUSSION 



We study a small black hole located on the TeV brane in Randall-Sundrum I (two branes) scenario. We expand the 
metric in e, which is the ratio between the Schwarzschild radius and the bulk curvature length. We find the solution 
up to first order in e. The solution satisfies Einstein's equations in the bulk and Israel junction conditions on the TeV 
brane. The asymptotic form of the metric is fixed by the weak field approximation (linearized gravity). 



In Randall-Sundrum II (single brane scenario) 6] it was conjectured that there were no large static black holes 
localized on the brane |l7j. However, small black holes might exist and one can apply our method of e expansion to 
find them. Although the linearized solution for a single brane is very different from that for two branes, in first order 
of e they are similar. The only difference is e — > — e. So, the results for the e solution and the discussion of regularity 
are equivalent in the two scenarios. 

A crucial issue is the post linearized form of the asymptotic solution. For Four dimensional black holes the 
post Newtonian behavior is given in terms of an expansion in integer powers of G^M/r (the only dimensionless 
combination). The equivalent in RS models would be an expansion in integer powers of G§M/(r£). However, if we 
assume that the post linearized metric includes only integer powers of the mass then the surface gravity on the horizon 
is not constant. This means that the horizon is singular, and the configuration of a small mass on the TeV brane 
describes a naked singularity. As we mentioned earlier, the singularity is naked from the bulk but it is covered on the 
brane. 

In section we have shown that there is a solution with a regular horizon, i.e. the surface gravity is constant and 
the Kretchmann scalar is finite up to first order in e. However, this requires a different asymptotic form of the metric, 
as one must include half integer powers of the mass. In general, the expansion parameter for the asymptotic solution is 
(G5M) 1 / 2 , but still the lowest order is (G5M) 1 and not (G5M) 1 / 2 otherwise there will not be a well defined conserved 
mass. Surprisingly, the first law of black hole thermodynamics, dM — T dS, eliminates the term of order (G5M) 3 / 2 , 
as well, from the asymptotic expansion. The leading term in the post Newtonian potential remains (G5M) 2 /(r 2 £ 2 ), 
which is the term responsible for the precession of perihelion calculations in the Schwarzschild metric. The next order 
term is (G5M) 5 / 2 . This term does not exist in the four dimensional Schwarzschild metric, but it is very hard to detect 
such a term in any asymptotic measurement. 

As we mentioned earlier, these results are also valid for small black holes in RS single brane scenario. In principle, 
the post Newtonian behavior can be verified numerically. However, the asymptotic behavior of the numerical solution 
of is not well known |18| . Future numerical studies with a regular horizon might detect the fractional power 
behavior in the post Newtonian potential. 

The configuration discussed in this paper, of a small (few TeV) mass on the brane in RS scenario, can be the final 
state of high energy particle collisions. In some sense such a collision is similar to the collapse of matter in ordinary 
four dimensional gravity, it might end up as a black hole or generate a naked singularity If the mass of the object 
is somewhat bigger than the Planck mass (TeV) the behavior of these two possibilities is very different. The black 
hole will radiate thermally. It will mainly radiate ordinary particles on the brane. The naked singularity has no 
temperature, it evaporates in an explosive fashion in the order of Planck time. It will mainly radiate gravitons into 
the bulk (where the singularity is naked) and not on the brane (where it is covered), further investigations of this 
issue will be carried out elsewhere. 



This work is supported in part by the U.S. Department of Energy Grant No. DE-FG02-84ER40153. We thank 
Richard Gass and Cenalo Vaz for fruitful discussions. 
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